A new Tolman test of a cosmic distance duality relation at 21 cm 
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Under certain general conditions in an expanding universe, tiie luminosity distance (di) and 
angular diameter distance (^a) are connected by the Etherington relation as d_t = dA{l + zf" . The 
Tolman test suggests the use of objects of known surface brightness, to test this relation. In this 
letter, we propose the use of redshifted 21 cm signal from disk galaxies, where neutral hydrogen (HI) 
masses are seen to be almost linearly correlated with surface area, to conduct a new Tolman test. 
We construct simulated catalogs of galaxies, with the observed size-luminosity relation and realistic 
redshift evolution of HI mass functions, likely to be detected with the planned Square Kilometer 
Array (SKA). We demonstrate that these observations may soon provide the best implementation 
of the Tolman test to detect any violation of the cosmic distance duality relation. 
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Introduction. — Observed fluxes and apparent sizes of 
cosmologically distant objects determine their dL and dA 
respectively, which are tied by the Etherington [l[ rela- 
tion asdi = dA{i + z)'^- This relation holds true for 
photons traveling on unique null geodesies in an expand- 
ing universe described by a metric theory of gravity and is 
guaranteed by phase space conservation of photons along 
with Lorentz invariance. Testing the Etherington rela- 
tion is important because any violation of it would be a 
smoking gun for new physics. At the same time, it is im- 
portant to distinguish between violations occurring due 
to new physics and astrophysical effects; often the signa- 
tures of these effects may not be distinguishable. Among 
the important astrophysical effects are dimming due to 
the intergalactic dust and gravitational lensing effects 2] ; 
while signatures of new physics could arise for example 
from interactions of photons with the dark sector. 

Tolman Q suggested that measurements of angular 
extension, luminosity and redshift of objects with a stan- 
dard surface brightness may be used to test this cosmic 
distance duality relation (CDDR) [28^1. In this letter, we 
propose 21 cm observations of disk galaxies as a new 
Tolman test that promises a systematics-free and clean 
probe of the CDDR. The advantage of our test over many 
of the existing ones is that it is immune to many astro- 
physical uncertainties. We show that future observations 
like those with the planned SKA would easily be able to 
probe this relation up to z ~ 1. We also discuss possi- 
ble strategies to obtain significantly tighter constraints 
on any possible violation of this relation. 

Existing Tests. — Sandage and Lubin [3] performed op- 
tical photometry of nearby galaxies and compared them 
with a sample of early- type galaxies at high redshift 0, 
performing the Tolman surface brightness test. The sur- 
face brightness was seen to be falling at a rate signifi- 
cantly slower than the oc (1 -I- z)"'' expected from the 
Etherington relation. However, the authors suggested 
that this discrepancy may be reconciled with an expand- 
ing universe by invoking luminosity evolution of early- 
type galaxies. 



Another way to test the CDDR is to make indepen- 
dent measurements of the cosmic distances and dA 
using different objects. References 0,0,01 provide upper 
limits on the cosmic transparency by observing the flux 
dimming of standard candles Q along with independent 
measurements of a cosmic distance scale through baryon 
acoustic oscillation measurements [9|| and H{z) [l^ data. 
However, all these methods suffer from uncertainties due 
to absorption of optical photons by dust; and also the 
possibility of being affected by lensing magnification ef- 
fects. 

Holanda et al. 



ll| have also examined the consistency 



of the CDDR through combined X-ray and SZ observa- 
tions of galaxy clusters. Although this method has the 
advantage of being unaffected by lensing magnification 
bias, their results depend largely on the cluster model 
and could suffer from unknown systematics. 

Our proposed test at 21 cm. — The Arecibo Dual-Beam 
burvey (ADBS) covered ~ 430 deg^ of the sky using the 
Arecibo main beam with a velocity coverage of —654 
to 7977 km s'^ [13] to detect the emission at 21 cm 
from galaxies. Out of the 265 identified galaxies, only 
a few were resolved in the Arecibo observations. Rosen- 
berg and Schneider [l3| made spatially resolved D-array 
follow-up observations with the Very Large Array (VLA) 
for 84 of these galaxies. Of these, 50 well resolved sources 
were used [ill to determine their HI masses and surface 
areas. In this study we use this sample to calibrate the 
size luminosity correlations to be used in the generation 
of our mock catalogs. 

HI masses and physical sizes of ADBS [l3| galaxies are 
found to be consistent with a nearly constant average HI 
surface density of the order ~ lO"^ kpc~^ [13]. The 
HI masses of individual ADBS galaxies, spanning more 
than three orders of magnitude, deviate by only ~ 0.13 
dex (Icr) from those expected from a constant HI surface 
density. This surprising relation has recently been ex- 
plained by Chakraborti 14] as the result of self-regulated 
star formation, driven by the competition between gravi- 
tational instabilities in a rotationally supported disk and 
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FIG. 1: The figure on the left represents the expected redshift distribution from some upcoming surveys. All angles in this 
figure are scaled up by a factor of 4 for representative purposes. The survey specifications and the la marginalized constraints 
on the parameters B, a and e are listed in the table on the right. 



mechanical feedback from supernovae. Silk fisj shows 
that supernovae drive up the porosity of the gas disk and 
stabilize it by increasing the gas dispersion velocity. Ana- 
lytic [in and simulation [3] results suggest that the gas 
dispersion velocity driven by supernova feedback is in- 
sensitive to the star formation rate (SFR) [2^ as well as 
metallicity [s^ . We find the redshift dependent effects in 
the evolution of the luminosity-area relation to be much 
lower than the intrinsic scatter in the HI masses. 

The present galaxy sample lies on the fundamental 
line |l4| even though the sample probes a large range 
in size and metallicity from dwarf galaxies to giant spi- 
rals. These objects with nearly constant surface bright- 
ness which is both documented 13] and understood .14], 
then provide us with astrophysical sources to give a di- 
rect measurement of the distance ratio dL/d-A at various 
redshifts. 

The HI emission starts out at 21 cm radio in the rest 
frame of the source. It gets redshifted out of the reso- 
nance line to longer wavelengths as the universe expands 
during its travel to the observer. The universe is essen- 
tially transparent at these wavelengths as their propaga- 
tion is not affected by gas and dust. Only HI can reabsorb 
the 21 cm emission before it is redshifted out of the res- 
onance line. But the absence of the Gunn and Peterson 
[r?! absorption feature for quasars below z ^ 6, confirm 
that there is hardly any HI (nni ^ 6 x 10^^^ atoms cm^"^) 
along these lines of sight. Considering the cross-section of 
the 21 cm hyperfine transition, we get the optical depth 
as ^ 10"^, which is negligible when compared to the in- 
trinsic scatter. This is a clear advantage over other tests 
of the CDDR at visible wavelengths where the effect of 
dimming due to dust is inseparable from effects due to 
new physics. Additionally, the use of the same objects to 
measure the ratio /dyi also means that there is no risk 
of lensing magnification bias, since the surface brightness 
of a single source remains preserved after gravitational 
lensing |18l |. 

Simulated Galaxy Catalogs. — To asses the potential of 
future radio observations in probing the CDDR using our 
proposed Tolman test, we create mock galaxy catalogs 



with a realistic redshift distribution that would be ex- 
pected to be seen in these surveys (see Fig. [ij. We 
briefly outline the procedure involved in doing so. 

Properties of the Galaxies: The luminosity, L and area, 
A of disk galaxies are seen to be related almost linearly 



L = 10^ {A/Aof Jy km s'^ Mpc^ 



(1) 



where we choose = 1 Mpc^. Fig. [5] shows the con- 
straints on the parameters B and a with the existing 
data of the ADBS galaxies. 

We use the Schechter mass function from Zwaan et al. 
[lit to model the HI mass distribution at low redshifts. 



dN 



r /Mhi 



dVdMm Af*jVAf*j 



/3 



exp 



Mhi 



HI 



(2) 



with /3 = -1.37 and log(Mj^i/M0) = 9.80. We hold 
Afj^j constant as it is likely to be regulated by the mass 
scale beyond which AGN feedback dominates supernova 
feedback f20j. However, as HI mass is directly related 
to the star formation rate, we choose a redshift depen- 
dent normalization for the HI mass function; 0*{z) = 
6 X 10-3(1 + 7.6z)/(l -I- (z/3.3)5-3), such that it matches 
the shape of SFR(z) from Hopkins and Beacom [21] . Any 
uncertainty in modeling the mass function here, only af- 
fects the number of galaxies yielded by the simulated 
survey and has no effect on the actual Tolman test. 

Following Rohlfs and Wilson [23], we relate the HI 
masses_to the flux / as, Mhi/Mq = 2.356 x 10^ /d|(l + 



z) ^ |3l|, where / is in Jy km s ^ and d^ is in Mpc. 



Corresponding to the the flux cut of a given survey, we 
calculate the limiting luminosities L(z) and correspond- 
ing masses M^\'^\z). The flux from the galaxy is as- 
sumed to be spread out, due to rotation, over a velocity 
range of Ay = 200 km s"^. If one assumes the weak 
scaling (AV oc M^^) from Briggs and Rao [l^] and in- 
clination effects, one would expect ~ 30% more sources 
[24I . In addition to the flux cut-off we also impose a size 
cut-off according to the beam size of the survey. The 
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FIG. 2: la & 2-cr marginalized constraints from 21 cm galaxy surveys - (Left) The green regions indicate parameter constraints 
obtained from observed ADBS galaxies (in which we fix e=0), while the orange regions represents constraints (marginalized 
over e) from a simulated galaxy catalog expected from the 5xeVLA survey. (Right) The green/blue regions correspond to the 
expected constraints from simulated galaxies in the SKA deep/wide surveys respectively. Note, how the degeneracy is broken 
(magenta) by combining a deep and wide survey. 



number of galaxies above the limiting mass of the survey 
with a coverage area Afi in a redshift bin Az is given by, 

dV f°° dN 
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(3) 



where V is the comoving volume upto a redshift z. This 
gives the luminosity and redshift distribution of galaxies 
and is used to create the simulated catalogs of galaxy 
data containing z, L, S and cr(S'), where a{S) obtains 
contribution mostly from the intrinsic scatter in the data 
as calibrated from the ADBS galaxies. 

Surveys at 21 cm: The properties of the simulated 
surveys are presented in Fig. [T] "ADBS + VLA" rep- 
resents the constraints on B and a obtained using the 
real data [l3| • "5 x eVLA" represents results from a sim- 
ulated catalog, assuming 80 hrs of on-source time using 
a hypothetical telescope with 5 times the sensitivity of 
the present Extended VLA (eVLA). This may soon be- 
come a reality with upgrades at the eVLA or upcoming 
pathfinder missions of the SKA. "SKA deep" represents 
a single snapshot (2 minutes) with the fully operational 
SKA. "SKA wide" represents the same amount of SKA 
time, but spent on 100 different fields instead of 1. While 
the SKA deep is good for constraining any redshift evolu- 
tion in the CDDR, the SKA wide is good for constraining 
the actual size luminosity relation of the galaxies. 

Results. — We modify the Etherington relation to test 
for any violation by using a simple parametrized form, 
dL/dA = (1 + Using the definitions d\iz) = A/S 

and c?|(z) = (1 + z)i/(47r/) [31], we relate the observed 
angular area (5) to observed flux (/) as follows. 



S{z,f) = 



47r/(l + z) 



3+2e \ - 



(4) 



,10Sd^(z)2(o- 

The almost linear relation (a ~ 1) between HI masses 
and surface area of disk galaxies implies that the residual 



cosmological dependence (through dA{z)) in the above 
equation is only expected to be weak. For the cosmology 
we assume the standard ACDM model fs^l keeping priors 
(from WMAP7 results [l^l) on ft. and We perform 
a MCMC likelihood analysis of the simulated data with 
the parameters - -B, a, e, h and flm- The likelihoods 
from the data are computed as exp(— x^/2), where = 



(log(5(z.,/.))-log(5,))M2 



We first analyze the ADBS data keeping the parameter 
e to be fixed, since the existing data is at very low red- 
shifts (z < 0.025) which is clearly insufficient to test for 
violations of the CDDR. We find the best fit parameters 
values as 5=9.003 and a=1.128 (see Fig. ^. We use 
these values for B and a along with the fiducial value of 
e~Q to create the mock catalog following the procedure 
as explained previously. 

We then analyze the simulated catalogs for some of 
the upcoming 21 cm observations like SxeVLA and the 
SKA. All of these surveys would throw up large number 
of galaxies in detection, which means that we may sim- 
ply self-calibrate the scaling relation of 21 cm luminosity 
with area (see eq. ([Ij). We discuss the prospects of the 
surveys to probe the CDDR by obtaining marginalized 
constraints on e, the parameter for violation of this rela- 
tion. 

The marginalized constrains on e from 5 x eVL A would 
be comparable to some of the present constraints [26| . 
with Ae=0.095. These constraints can be significantly 
improved by the next generation telescopes like the SKA. 
The SKA deep survey with a very good flux sensitivity is 
expected to detect galaxies upto a high redshift of z ^ 1. 
The 100 sq. deg. wide SKA survey would have many 
more objects, however, because of its higher flux cut, it 
would mostly detect low mass galaxies at lower redshifts, 
z ~ 0.4 (see Fig. [T]). Both the deep and wide surveys 
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from the SKA would give Ae=0.026; however, because 
these surveys probe the relation in eq. ([Ij at very differ- 
ent range of masses they place constraints along different 
directions. This fact can be used to a great advantage 
by using the deep and the wide surveys in conjunction; 
as shown in Fig. [5J this is excellent at breaking the de- 
generacy between e and a. This wedding cake 27[ survey 
strategy would strongly constrain e to 0.01 while also 
constraining the relation in eq. ([T]) to a high precision 
with AB=0.007 and Aa=0.002. 

Discussions. — The planned SKA will allow the detec- 
tion of 21 cm emission from distant {z ~ 1) galaxies with 
small observing times. Apart from being useful in per- 
forming the Tolman test, the same data will be useful 
for cosmology and galaxy evolution studies. The deep- 
est observations may allow one to reach z ^ 6 and start 
probing the HI content of the first galaxies, which would 
have played a vital role in the epoch of reionization. 

Self regulation of mechanical feedback from supernovae 
and gravitational instability of the star forming disks 
keeps the surface density nearly constant, irrespective of 
the global star formation rates. The almost constant ob- 
served surface density of HI in disk galaxies makes them 
ideal sources for the Tolman test - Firstly, by using the 
same objects to directly probe the ratio d^/dA at vari- 
ous redshifts, the results are not affected by many possi- 
ble astrophysical model uncertainties. Secondly, our test 
also promises to be immune to bias from effects of grav- 
itational lensing as it preserves the surface brightness of 
such sources. Lastly, our proposed test does not fix the 
scaling relation between luminosity and the area of HI 
regions in galaxies; the upcoming surveys would simply 
detect enough number of these galaxies to permit a self- 
calibration of this relation. 

We have shown that observations at the 21 cm wave- 
length opens up the exciting possibility of testing the 
CDDR in a completely independent manner. Probing 
this relation with our proposed method, offers several 
advantages over all of the existing tests. At radio fre- 
quencies the universe is essentially transparent. Thus 
one can cleanly distinguish between violations of this re- 
lation arising from astrophysical effects and those that 
may arise from new physics. Our proposed test of the 
CDDR may soon enable us to probe the possibility of 
photons coupling to the dark sector. 
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